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The hierarchy of exact equations is given that relates two-spatial-point velocity structure func- 
tions of arbitrary order with other statistics. Because no assumption is used, the exact statistical 
equations can apply to any flow for which the Navier-Stokes equations are accurate, and they apply 
no matter how small the number of samples in the ensemble. The exact statistical equations can 
be used to verify DNS computations and to detect their limitations. For example, if DNS data are 
used to evaluate the exact statistical equations, then the equations should balance to within numer- 
ical precision, otherwise a computational problem is indicated. The equations allow quantification 
of the approach to local homogeneity and to local isotropy. Testing the balance of the equations 
allows detection of scaling ranges for quantification of scaling-range exponents. The second-order 
equations lead to Kolmogorov's equation. All higher-order equations contain a statistic composed 
of one factor of the two-point difference of the pressure gradient multiplied by factors of velocity 
difference. Investigation of this pressure-gradient-difference statistic can reveal much about two 
issues: 1) whether or not different components of the velocity structure function of given order 
have differing exponents in the inertial range, and 2) the increasing deviation of those exponents 
from Kolmogorov scaling as the order inc reases. Full disclosure of the mathematical methods is in 
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I. INTRODUCTION 



Kolmogorov's (1941) equation and Yaglom's equation were the first two equations of the "dynamic theory" 
of the local structure of turbulence. The name "dynamic theory" was originated by Monin & Yaglom (1975) (their 
Sec. 22) to mean the derivation of equations relating structure functions by use of the Navier-Stokes equation and/or 
the scalar conservation equation, and the investigation of the resulting statistical equations. Monin & Yaglom 
(1975) pointed out that the dynamic theory gives important relationships among structure functions, and that these 
relationships provide extensions of predictions based on dimensional analysis. Theoretical studies (Lindborg 1996; 
Hill 1997a) clarified the assumptions that are the basis of Kolmogorov's equation and give equations that are valid 
for anisotropic and locally homogeneous turbulence as well as for the case of local isotropy and local homogeneity. 
Antonia et al. (1983) and Chambers & Antonia (1984) used experimental data to study of the balance of the classic 
equations of Kolmogorov and Yaglom. There is renewed interest in examining the balance of those equations using 
both experimental and DNS data and in generalizing the equations to cases of inhomogeneous, nonstationary, and 
anisotropic turbulence (Lindborg 1999; Danaila et al, 1999a, b,c; Antonia et al. 2000). Whereas Kolmogorov's 
(1941) equation relates 2nd- and 3rd-order velocity structure functions, the next-order dynamic equation relates 3rd- 
and 4th-order structure functions and a pressure-gradient, velocity- velocity structure function. The balance of that 
next-order equation has been examined by means of experimental and DNS data; this showed the behavior of the 
pressure-gradient, velocity-velocity structure function (Hill & Boratav 2001). There is now interest in dynamic-theory 
equations of arbitrarily high order A" (Yakhot 2001). Such equations relate velocity structure functions of order N 
and AT + 1 and other statistics. Those equations are given in this paper. 

Using the assumptions of local homogeneity, local isotropy and the Navier-Stokes equation, Yakhot (2001) 
derived the equation for the characteristic function of the probability distribution of two-point velocity differences. 
He uses that equation to derive higher-order dynamic equations. Equations for arbitrarily high-order structure 
functions can be obtained by repeated application of his differentiation procedure. Yakhot (2001) studies the inertial- 
range, deduces a closure, and thereby determines the incrtial-range scaling exponents of velocity structure functions. 
Yakhot's study is the first to make significant use of dynamic-theory equations to determine scaling exponents. 

The purposes and theoretical method of the present paper differ from those of Yakhot (2001), but one purpose 
is to verify Yakhot's equations from our distinctly different derivation. That verification is given in Sec. 5. In 
Sec. 2, exact statistical equations relating velocity structure functions of any order are derived from the Navier- 
Stokes equation. "Exact" means that no assumptions are made other than the assumption that the Navier-Stokes 
equation and incompressibility are accurate. Since the equations are exact, they apply to any flow, including laminar 
flow and inhomogeneous and anisotropic turbulent flow. The exact statistical equations can be used to verify DNS 
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computations and detect their limitations. New experimental methods of Dahm and colleagues (Su & Dahm 1996) 
can also be tested. For example, if DNS data are used to evaluate the exact statistical equations, then the equations 
should balance to within numerical precision, otherwise a computational problem is indicated. In Sec. 3, statistical 
equations valid for locally homogeneous and anisotropic turbulence are obtained from the exact equations; those 
equations can be used with DNS or experimental (Su & Dahm 1996) data to study the approach to local homogeneity 
of a particular flow. This can be done by quantifying the terms that are neglected when passing from exact equations 
to the locally homogeneous case, and by quantifying changes in the retained terms as local homogeneity is approached 
when the spatial separation vector is decreased. In Sec. 4, statistical equations valid for locally isotropic and locally 
homogeneous turbulence are obtained from those for the locally homogeneous case. The approach to local isotropy 
can be studied by means analogous to the above described evaluation of local homogeneity. Such studies might 
shed light on the observed persistence of anisotropy (Pumir & Shraiman 1995; Shen & Warhaft 2000). All dynamic- 
theory equations are now available to extend the above-mentioned previous studies of the balance of dynamic-theory 
equations. 

Incompressibility requires that the different components of the second-order velocity structure function have 
the same scaling exponent in the inertial range. The same is true for the third-order structure function. However, 
at fourth and higher order there is no such requirement. There have been many studies of the possibility that the 
inertial-range scaling exponents of the various structure-function components are unequal (e.g., Chen et al. 1997; 
Boratav & Pelz 1997; Boratav 1997; Grossmann et al. 1997; van de Water & Herweijer 1999; Camussi & Benzi 1997; 
Dhruva et al 1997; Antonia et al. 1998; Kahaleras et al. 1996; Noullez et al. 1997; Nelkin 1999; Zhou & Antonia 
2000; Kerr et al. 2001). The usefulness of applying the higher-order dynamic-theory equations to those investigations 
is considered in Sec. 6. 

Derivation of the equations produces substantial mathematical detail. Matrix-based algorithms are invented 
such that the isotropic formulas for the divergence and Laplacian of isotropic tensors of any order can be generated 
by computer. The details of this mathematics are available and are herein referred to as the Archive.^ 

II. EXACT TWO-POINT EQUATIONS 

The Navier-Stokes equation for velocity component Ui (x, t) is 

d t Ui(x,t) +m„(x, t)d Xn Ui{x, t) = -d Xi p(x,t) + vd Xn d Xn Ui(x, t), (1) 

and the incompressibility condition is d Xn u n (x.,t) =0. In (|l|), p(x, t) is the pressure divided by the density (density 
is constant), v is kinematic viscosity, and d denotes partial differentiation with respect to its subscript variable. 
Summation is implied by repeated Roman indices. Consider another point x' such that x' and x are independent 
variables. For brevity, let Ui = «i(x, t), u[ = iti(x',i), etc. Require that x and x' have no relative motion. Then 
d Xi u'j = 0, d x 'Uj = 0, etc., and d t is performed with both x and x' fixed. The change of independent variables from 
x and x' to the sum and difference independent variables is: 

X = (x + x') /2 and r = x x', and define r = |r| . (2) 

The relationship between the partial derivatives is 

d Xl = d n + \dxi , d x >. = -d ri + ^d Xl , d Xi = d Xi + d x ^ , d r% = ^ (d Xi - d x ^j . (3) 

The change of variables organizes the equations in a revealing way because of the following properties. In the case of 
homogeneous turbulence, dxi operating on a statistic produces zero because that derivative is the rate of change with 
respect to the place where the measurement is performed. Consider a term in an equation composed of dx t operating 
on a statistic. For locally homogeneous turbulence, that term becomes negligible as r is decreased relative to the 
integral scale. For the homogeneous and locally homogeneous cases, the statistical equations retain their dependence 
on r, which is the displacement vector of two points of measurement. Subtracting ([!]) at x' from ([!]) at x, performing 
the change of variables (^|), and using (^) gives 
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d t Vi + U n d Xn v l + v n d rn v l = -Pi + v (d Xn d Xn Vi + d x > n d x / n Vi) , 

where v l = u i - u-, U n = (w„ + u' n ) /2, Pi = d Xi p - d x rp . 



(4) 
(5) 



Now multiply (|J) by the product VjV k ■ ■ ■ vi which contains N — 1 factors of velocity difference, each factor 
having a distinct index. Sum the N equations as required to produce symmetry under interchange of each pair of 
indices, excluding the summation index n. French braces, i.e., {o}, denote the sum of all terms of a given type that 
produce symmetry under interchange of each pair of indices. The differentiation chain rule gives 



{vjv k ■ ■ ■ vid t Vi} = d t {vjV k ■ ■ ■ viVi) , 
{VjVk ■ ■ ■ viU n d x „Vi} = U n d x „ (vjv k ■ ■ ■ vm) = d Xn [U n VjVk ■ ■ ■ vm 

{VjV k ■ ■ ■ VlV n d rn Vi} = V n d Tn (VjV k ■ ■ • = 9 r „ {v n V.jV k ■ ■ ■ VlVi) . 



(6) 
(7) 
(8) 



The right-most expressions in (Q) and (0) follow from the incompressibility property obtained from (j|) and the fact 
that d Xi u'j = 0, d x >iij = 0; namely, dx„U n = 0,dx n v n — 0,d rn U n = 0,d rn v n = 0. The viscous term in (^) produces 
v {vjV k ■ ■ ■ vi (d Xn d Xrl Vi + d x i d x i Vi) }; this expression is treated in the Archive. Thereby 



d t (vj ■ ■ ■ Vi) + d Xn (U n Vj ■••«*) + dr n (v n Vj ■ ■ ■ vi) 

1 



{vj ■ ■ ■ viPi} + 2v 



d rrh dr n + -d Xn d Xn (vj ■ ■ ■ - {v k ■ ■ ■ vieij} 



(9) 



where e i3 = (d Xn Ui) (d Xn Uj) + (d x > n u[) (d x > n Uj) = (d Xn Vi) (d Xn Vj) + (d x > n Vi) {d x > n Vj 



(10) 



The quantity {vj ■ ■ • UjPj} can be expressed differently on the basis that allows to be written as Pi = 
dxi (p — p')- The derivation is in the Archive; the alternative formula is 



{vjv k ■ ■ ■ viPi] = {d Xl [vjV k ■ ■ ■ vi (p-p')]} -(N-l) ip-p') {(( 
where the rate of strain tensor Sij is defined by sy = [d Xi Uj + d Xj Ui) /2. 



vi) 



(11) 



A. Hierarchy of exact statistical equations 



Consider the ensemble average because it commutes with temporal and spatial derivatives. The above notation 
of explicit indices is burdensome. Because the tensors are symmetric, it suffices to show only the number of indices. 
Define the following statistical tensors which are symmetric under interchange of any pair of indices, excluding the 



summation index n in the definition of F 



D 



[N] 



( V 3 



"i) > F [JV+1] 



[N+l]- 

(U n Vj 



Vi) , T[AT] 



({vj ■ ■ ■ viPi}) , E[jvj = ({v k ■ ■ ■ vi^}) 



(12) 



where angle brackets () denote the ensemble average, and the subscripts N and N + l within square brackets denote 
the number of indices. The left-hand side of each definition in ( |l2|) is in implicit-index notation for which only the 
number of indices is given; the right-hand sides in (|l^) are in explicit-index notation. The argument list for each 
tensor is understood to be (X, r, t). The ensemble average of (0) is 



[N] 



Vx»F 



[N+l] 



V r «D 



[N+l] 



■[N] 



2v 



-V 



x 



D 



[N] 



E 



[N] 



(13) 



where, Vx • F[jv+i] = dx n (U n Vj ■ ■ ■ v%) , V r • D[ Ar+1 ] = d rn {v n Vj •••«,•}, = d rn d rn , V x = dx n dx n ■ The notations 
Vx», Vxj V r », and are the divergence and Laplacian operators in X-space and r-space, respectively. 



III. HOMOGENEOUS AND LOCALLY HOMOGENEOUS TURBULENCE 



Consider homogeneous turbulence and locally homogeneous turbulence; the latter applies for small r and large 
Reynolds number. The variation of the statistics with the location of measurement or of evaluation is zero for the 
homogeneous case and is neglected for the locally homogeneous case. Since that location is X, the result of Vx 
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operating on a statistic vanishes or is neglected as the case may be. Thus the terms Vx • F[jv+i] and |V x D[jv] are 
deleted in (|l3|); then ( |l3| ) becomes, 

d t B [N] + V r • D [JV+1] = -T [N] + 2v [V*D [JV] - E[jv]] . (14) 

Because the X-dependence is deleted, the argument list for each tensor is understood to be (r, t). Note 
that dtDm] is n °t necessarily negligible for homogeneous turbulence. The ensemble average of (O) contains 
{{dxt [vjVk ■ ■ ■ vi (p — p')]}) = {dxi (vjVk ■ ■ ■ vi(p — p')}} = {0} = 0. Thus, ( |TT| ) gives the alternative that 

T im = -(N-l)((p-p'){(s ij -s' ij )v k ---v l }). (15) 

One distinction between ( [l4| ) and the hierarchy equations given in equations (13) and (17) by Arad et al. (1999) is 
that their t- and r-derivatives operate on only one velocity difference within their product of such differences, whereas 
the derivatives in (0) and thus in (|l4|) operate on all N of the velocity differences. 



IV. ISOTROPIC AND LOCALLY ISOTROPIC TURBULENCE 



Consider isotropic turbulence and locally isotropic turbulence; the latter applies for small r and large Reynolds 
number. Locally isotropic flows are a subset of locally homogeneous flows (Monin & Yaglom Sec. 13.3, 1975) and 
similarly for the relationship between isotropic and homogeneous flows. Thus, the dynamical equations for locally 
isotropic and isotropic turbulence are obtained from (|l4|) such that the variable X and the terms Vx • F^+i] and 
^V x D[jv] (see |l3|) do not appear. The tensors D[jvj, T[n]i and E[jy] in ( O ) obey the isotropic formula. The 
Kronecker deltaJy is defined by <5,y = 1 if i = j and <5y = if i ^ j. Let dnp] denote the product of P Kronecker 
deltas having 2P distinct indices, and let WW] (r) denote the product of N factors a each with a distinct index; 
the argument r is omitted when clarity does not suffer. Because each tensor in ( |l2] ) is symmetric under interchange 
of any two indices, their isotropic formulas are particularly simple. Each formula is a sum of M + 1 terms where 
M = N/2 if N is even, and M = (N — 1) /2 if N is odd. Each term is the product of a distinct scalar function with 
a Wjjv] and a 5pp]. From one term to the next a pair of indices is transferred from a WW] to a 5[2P\) examples are 
in the Archive. For the tensor D^j, denote the Pth scalar function by Dn,p (r, t). The isotropic formula for T>[n] is 



T> [N] (r,i) =J2 Dn,p (r,t) {W [N _ 2P] (r) 6 [2P] } , (16) 



p=o 



and the isotropic formulas for Trjy] and Erjy] have the analogous notation. Recall from Sec. 2 the meaning of the 
notation {o} whereby {W[ N _ 2 p] (r) fcp]} denotes the sum of all terms of the type W[ A r_ 2 p] ( r ) <5[2P] that produce 
symmetry under interchange of each pair of indices. An example is {W^j (r) <5[2]} = + -jr&H + ^$jk- 

A special Cartesian coordinate system simplifies the isotropic formulas. This coordinate system has the positive 
1-axis parallel to the direction of r, and the 2- and 3- axes are therefore perpendicular to r. Let N%, N 2 , and N$ be the 
number of indices of a component of Dr/v] that are 1,2, and 3, respectively; such that N = N\ + N 2 + N% . Because of 
symmetry, the order of indices is immaterial so that a component of Dm can be identified by Ni, N 2 , and N%. Thus, 
denote a component of T)[n] by ^[jVi.JVa.iVs] which is a function of r and t. Likewise, { W[jy-2P] ( r ) <^[2P] } < Nl N2 N ^ 1S a 
specific component of the tensor {W[ A r_ 2 p] ( r ) <5[2P] }• If; m ( |l6|) N\ of the indices are assigned the value 1, and N 2 and 
-/V3 of the indices are assigned the values 2 and 3, respectively, then f [jVi.JVa.iVs] and {~W[n-2P] ( r ) ^[2P]}r Wl N2 N ^ wn ^ 
appear on the left-hand and right-hand sides of respectively. The {W[ A r_ 2 p] ( r ) $[2P]}[ Nl N ^ N i are numerical 
coefficients that do not depend on r because & = - = 1, ^ = ^ = 0. From the Archive, the values of the coefficients 
are: 

if 2P < N 2 + N 3 then {W [A r_ 2P] (r) S\2P\}[N lt N a ,N,] = °' otherwise ' ( 17 ) 

N 2 \,fN 3 \,f n ^2 A^3 



{W [Ar _ 2 pj (r) S m } [NuN2>N3] = N^miNsl/ [(JV - 2P)!2^ j ! j ! \p - tf - ^ j !j . (18) 

By applying ( |l6[]l8| ) for all combinations of indices, one can determine which components D[n 1 .n 2 ,n 3 \ are zero 
and which are nonzero, identify M + 1 linearly independent equations that determine the Dn p in terms of M + 1 of 
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the DtN lt N 2> N 3 ], and find algebraic relationships between the remaining nonzero D^ lt j^ 2t j^ 3 j. The derivations are in 
the Archive; a summary follows. 

A component -D[7Vi,Af 2 ,Af 3 ] is nonzero only if both N% and N$ are even, and therefore N\ is odd if N is odd, and 
Ni is even if N is even. Thereby, (M + 1) (M + 2) /2 components are nonzero. There are 3^ components of Drm; 
thus the other 3^ — (M + 1) [M + 2) /2 components are zero. 

There exist (M + l)M/2 kinematic relationships among the nonzero components of D[jyi. For each of the 
M + 1 cases of Ni , these relationships are expressed by the proportionality 



^[JVi,2L,0] : -D[A r i,2L-2,2] : £*[jVi,2£-4,4] ^[JVi,0, 2L] — 

[(2L)!0!/L!0!] : [{2L - 2)!2!/ (L - l)!l!] : [(2L - 4)!4!/ (L - 2)!2!] 



[0! (2L)!/0!i! 



(19) 



For N = 4 with L = 2, ( |19[) gives £^[o,4,o] : ^[0,2,2] : -O[o,o,4] = 12 : 4 : 12. In explicit-index notation this can be written 
as D 2222 = 3Z?2233 = -D3333, which was discovered by Millionshtchikov (1941) and is the only previously known such 
relationship. Now, all such relationships are known from (|l9|). 

There remain M + 1 linearly independent nonzero components of Dry]. This must be so because there are 
M + 1 terms in (fig) and the M + 1 scalar functions Dn,p therein must be related to M + 1 components. Consider 



the M + 1 linearly independent equations that determine the p in terms of M + 1 of the Ow^aTj,^]' F° r 
simplicity, the chosen components can all have N3 — 0; i.e., the choice of linearly independent components can be 
D [jv,o,0]> ^[jv-2,2,0]; ^[iV-4,4,0] j ' ' ') D [N-2M,2M,o] ■ As described above, assigning index values in ( |T6| ) results in the 
chosen components on the left-hand side and algebraic expressions on the right-hand side that contain the coefficients 
(00)" ^ n ^ e Archive, those equations are expressed in matrix form and solved by matrix inversion methods. Given 
experimental or DNS data or a theoretical formula for the chosen components, the solution of the algebraic equations 
determines the functions Djy^p in jl^); then ( |l6|) completely specifies the tensor Dry]. 

The matrix algorithm in the Archive is an efficient means of determining isotropic expressions for the terms 



D 



[iV+l] 



and V^Djjv] in (|14|). From the example for N = 2 in the Archive, (|by) gives the two scalar equations 



d t D u 



Dm -D122 — —T\ 



2v 



2 i 4 
-d r - -5 



= 2v 



:D 2 



d 2 D lL + -d r D u + -j (D 22 - D n ) 



2vE\ 



4e/3, 



(20) 



d t D 22 + [d r + - D 122 = -T 22 + 2u 



2 



d l r + -d r - — } D 22 
r r A 



= 2v 



d 2 r D 22 + -d r D 22 - — {D 22 - D n ) 



- 2vE 22 
4e/3, 



(21) 



where use was made of the fact (Hill 1997a) that local isotropy gives T\\ = T 22 = and 2vEu — 2vE 22 — 4e/3, where 
e is the average energy dissipation rate. Since (^p|-pTl) are the same as equations (43-44) of Hill (1997a), and since 
Hill (1997a) shows how these equations lead to Kolmogorov's equation and his 4/5 law, further discussion of ( p0f|2l| ) 
is unnecessary. From the example for N = 3 in the Archive, 



d t D u 



<>, + - I Dun - -D 1122 = -Tm + 2v [(V?D) m - E m ] 



<h Dvii y \ r + - ) D1122 - ^2222 = -T 122 + 2v [(V?D) 122 - E 122 ] 



(22) 
(23) 



(V»D) m ee [d 2 r + \d r Dm + y 2 D 122 = (-1 + \d r + d^j Dm, (24) 

(V"D) 122 = Idiu + (d 2 r + ~9 r - D X22 = i (1 - U r + 5d 2 + rd*) Dm- (25) 

The incompressibility condition, D\ 22 — g (Dm + rd r Dm), was used in (P^]25| ). Since Hill & Boratav (2001) discuss 
these equations and evaluate them using data, further discussion of (p2]-[25P is unnecessary. 

The terms dtD[N], — T[at], and — 2i/E[jv] in ( p^ ) have a repetitive structure in the isotropic equations; e.g., for 
N = 4 the 3 equations are 
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d t £mi + (V r .D [5] ) im =-r ull +2^[(V^D [4] ) mi - J B lul ], (26) 
a t ^ii 22 + (V r .D [5] ) 1122 =-T U22 + 2^[(V^D [4] ) 1122 -£; U22 ] , (27) 
d t D 2222 + (V r . D [5] ) 2222 = -T 2222 + 2u [(V?D [4] ) 2222 - E 2222 ] . (28) 

Thus, it suffices to give the isotropic formulas for the divergence V r • D[jv + i] and Laplacian V^D^; for N = 4 to 7, 
those isotropic formulas are given in Table 1. For N = 4 and 5 there are M + 1 = 3 equations; there are M + 1 = 4 
equations for both N = 6 and 7. 

{pr + f:) #[5,0,0] - ;:#[3,2,0] (<% + f^r ~ pr) #[4,0,0] + pr#[2,2,0] + ^3#[0,4,0] 
[pr + § ) #[3,2,0] - 37#[1,4,0] ^-D^O.O] + ( - + <9, 2 + frcU #[2, 2 ,0] + §fl#[0,4,0] 



<9 r + J. ) #[1,4,0] fTT#[2,2,0] + (- 3^T + + r^r) #[0,4,0] 

(<9 r + ~) #[6,0,0] - ^T#[4,2,0] (pr + f ® r ~ w) ^[5,0,01 _ p#[3, 2,0] + ft#[l,4,0] 

(9,- + |) #[4,2,0] - 7-0(2,4,0] fTT#[5,0,0] + ( — + + f^r) #[3,2,0] + ^3#[1,4,0] 

(Or + I) #[2,4,0] - £#[0,6,0] 57^[3,2,0] + + 9? + ^r) £> [M)0] 

f« _i_ 2 \ n 12 r> /a2 , 2i 12\ n i08 n , 920 n 416 n 

(Cr + f ) -^[7,0,0] - —^[5,2,0] [O r + ~ dr ~ ~ ) #[6,0,0] ~ —^[4,2,0] + 372^12,4,0] - 1572 #[0,6,0] 

(9 r + I) #[5,2,0] - ||#[3,4,0] fTT#[6,0,0] + (-57! + ^ + f #[4,2,0] + ^3#[2,4,0] ~ fpT#[0,6,0] 

57#[1,6,0] 5^#[4,2,0] + (~^r + d l + l d r) #[2,4,0] + 3^ 



(d r + -) #[1,6,0] 377#[2,4,o] + (^TsT 5 ' + 9 r + -d r ) #[0,6,0] 

(dr + |) #[8,0,0] - ^r#[6,2,0] (<% + ~ ^\^ 7 '°^~ ^^[5, 2,0] + ^^#J3,4,0] ~ ^T#[l,6,0] 

(pr + g) #[6,2,0] " |g#[4,4,0] f|#[7,0,0] + ( _ 2^|+ + 7®r) #[5,2,0] + i^^I 3 ' 4 ' ! ~~ Ipr#[l,6,0] 
I + l\ D[4A ° ] ~ f D P fi ^ 7| 2 '#[5,2,0] + ("l^g* + + \ dr \ D l 3 ^°] + T7^#[l:6,0] 
\8r + -) #[2,6,0] - vv#[0,8,0] 7y2-#[3,4,0] + (-35^ + ®r + T^r) #[1,6,0] 

Table 1: The isotropic formulas for V r • Drjy+1] are on the left and those for V^D^] are on the right. Cases 
N = A and 5 are the top 3 and second 3 rows, respectively. Cases N = 6 and 7 are the third 4 and bottom 4 rows, 
respectively. 



V. COMPARISON WITH PREVIOUS RESULTS 

The expression (d r + -) Dm — ^D\ 22 in ( po| ) is the same as equation (9) of Yakhot (2001), and (41) of Hill 
(1997a). The expression (d r + |) #im — f#n 22 in ( J2^ ) is the same as in the equation that follows Yakhot's eq.ll , and 
in equation (16) of Hill & Boratav (2001) and in equation (8) of Kurien & Sreenivasan (2001); (d r + *) #1122 — £#2222 
in ( p3| ) is the same as in equation (13) of Hill & Boratav (2001) and equation (10) of Kurien & Sreenivasan (2001). 
The expressions (d r + |) #[6,0,0] ~ ^#[4,2,0] and (d r + |) #[ 2>4 ,o] - ^#[0,6,0] for the case N = 5 in Table 1 are 
the same as in equations (9) and (10) of Kurien & Sreenivasan (2001). More generally, the isotropic formulas for 

V r • Drjv+i] for the case Ni — N, N 2 — N 3 — are (<9 r + |) f[jv,o,o] — 2 S~ #[JV-2,2,o] which agrees with the 
left-hand side of equation (7) of Yakhot (2001). The other components of V r • D[jy + i] were not given by Yakhot 
(2001). The expressions from the Laplacian in (20-21]) are the same as in (41-42) of Hill (1997a); and ( pi| - |25| ) are the 



same as (7-8) of Hill & Boratav (2001). All of the remaining results do not appear to have been given previously. 
The above comparisons are sufficient to verify the matrix algorithm for generating the structure-function equations 
to any desired order, as well as to independently validate the derivation of Yakhot (2001). 



VI. SUMMARY AND DISCUSSION 



The third paragraph of the introduction summarizes part of this paper and is not repeated here. In addition: 
All of the kinematic relationships (|l9|) between components of isotropic, symmetric structure functions of arbitrary 
order have been identified, whereas previously only one was known. All of the components that are zero have been 
identified (a recent experimental evaluation of some of them is given by Kurien & Sreenivasan 2000). The kinematic 
relationships show that the scaling exponents of certain different components must be equal; if the exponents are not 
equal when evaluated using one's data, then the kinematic relationships ( |l9| ) provide a measure of either the error 
in the exponents or the deviation from local isotropy. The dynamic equations of order N can be used to test the 
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extent of a scaling range for evaluation of scaling exponents of velocity structure functions of order N + 1 because the 
time-derivative and viscous terms should be zero in an inertial range. The graphical presentations of the balance of 
Kolmogorov's equation by Antonia et al. (1983), Chambers & Antonia (1984), Danaila et al. (1999a, b), and Antonia 
et al. (2000) show the extent of, or deviation from, inertial-range exponents. The higher-order equations given here 
can be used in an analogous manner. 

The energy dissipation rate e plays an essential role at all r in Kolmogorov's equation. In our formulation e 
) from the tensor components 1vE\\ and IvEii- On the other hand, for the next-order equations ( p2[ - 
|23| ) Hill (1997b) showed that the corresponding terms 2vEm and 1vE\ii are negligible in the inertial range. Yakhot 
(2001) shows that the components E[ N ] are negligible in the inertial range for all of the higher-order equations for 
which N is odd. Kolmogorov's (1941) inertial-range scaling using e and r as the only relevant parameters can be 
used to estimate the relative magnitudes of the term V r • D^+i] in (|l4|) to the terms 2v\/^'D[ N ] and 2u~E,y N y Doing 
so, the ratio of any nonzero component of 2z^V^D[jy] or 2^E[jy] to the corresponding component of V r • D[tv+i] is 

proportional to ^/r 4 / 3 ^ 1 / 3 = (r/r/) 4 ^ 3 , which asymptotically vanishes in the inertial range (n = (^ 3 /e) 1 ^ 4 ). Thus, 
both terms 2^VpD[jv] and 2^E[jy] are to be neglected in an inertial range if N > 2. 

One concludes that all equations of order higher than Kolmogorov's equation reduce to the isotropic formula 
for V r • D[jv+i] = — T[tv] in the inertial range. This formula shows that Tnyi 1S a t the heart of two issues that 
have received much attention: 1) whether or not different components of the velocity structure function D[jy+i] have 
differing exponents in the inertial range, and 2) the increasing deviation of those exponents from Kolmogorov scaling 
as N increases. The physical basis for the importance of T\n] is the importance of vortex tubes to the intermittency 
phenomenon (Pullin & Saffman 1998) combined with the fact that the pressure-gradient force is essential to the 
existence of vortex tubes; the pressure-gradient force prevents a vortex tube from cavitating despite the centrifugal 
force. Pressure gradients are the sinews of vortex tubes. Direct investigation of T[n] using DNS can reveal much 
about the two issues. 
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